Confinement Phase in Carbon-Nanotubes and the Extended Massive Schwinger Model 
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Carbon nanotube with electric fluxes confined in one dimension is studied. We show that a Coulomb 
interaction oc \x\ leads to a confinement phase with many properties similar to QCD in 4D. Low- 
energy physics is described by the massive Schwinger model with multi-species fermions labeled 
by the band and valley indices. We propose two means to detect this state. One is through an 
optical measurement of the exciton spectrum, which has been calculated via the 't Hooft-Berknoff 
equation with the light-front field theory. We show that the Gell-Mann— Oakes— Renner relation is 
satisfied by a dark exciton. The second is the nonlinear transport which is related to Coleman's 
"half-asymptotic" state. 
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Even after two decades from its discovery, fascination 
continues with the carbon nanotube (l|. One line of re- 
search is motivated by its possible application as opto- 
electronic devices. Another is more academic, in which 
possibilities are explored to realize new states of matter. 
In low dimensions quantum fluctuations are enhanced, 
which makes the nanotube an ideal host for strongly cor- 
related and clean quantum systems. One monumental 
result along this line is the Tomonaga-Luttinger liquid 
in single- wall metallic nanotubes [2H4j. I n the present 
work, we theoretically predict a realization of yet an- 
other interesting state, namely the confinement phase, in 
nanotubes. In this state, charged particles cannot exist 
as free asymptotic states where all the excitations are 
neutral bound states (i.e., excitons in the present case). 

A most famous realization of the confinement phase 
is the hadronic system where the effective quantum field 
theory is QCD in (3+l)D. By contrast, a carbon nan- 
otube is a one-dimensional object, and no dynamical non- 
abelian gauge field exists. Nonetheless, we propose here 
a confinement phase may exist; A key is the dimension- 
ality. In a seminal paper by Schwinger, a (1+1)D ver- 
sion of QED, i.e. the Schwinger model, was studiedQ. 
In this model, Dirac particles interact through a one- 
dimensional Coulomb potential (ex |x|). The long-range 
interaction mimics the confinement potential, and the 
groundstate is indeed in a confinement phase. The 
Schwinger model and its extensions were studied as a toy 
model of (3+l)D QCD to understand non-perturbative 
aspects of the confinement phase |6l-[l3|. However, up 
to now, no physical realization of this model is known. 
Our claim here is that the low-energy effective model 
for electrons in a nanotube can be, in certain situations, 
modeled by an extended massive Schwinger model. It 
is well known that a multi-species Dirac spectrum is re- 
alized in the band structure of nanotube within the ef- 
fective mass approximation. In order to realize the one- 
dimensional Coulomb potential, the electric fluxes must 
be confined along the tube. The situation with shielded 
ID electric flux can possibly be realized in (a) multi- 
wall nanotubes with metallic outer tubes, (b) a tube sur- 
rounded by metallic tubes, or (c) tubes embedded in a 



superconductor. While one may wonder if such shielding 
can really take place, there is a prominent example in 
two-dimensional (2D) organic crystals studied by Yam- 
aguchi et aZ.JlJ]. They have experimentally shown that 
a layered structure with a large difference in the dielec- 
tric constants leads to confinement of electric fluxes in 
the 2D plane where the Coulomb potential becomes log- 
arithmic. In this experiment, the long-range Coulomb 
interaction leads to a power- law current-electric field (J- 
E) characteristics with a temperature dependent power. 
This gives us a strong motivation to study nanotube with 
a ID Coulomb potential. The properties of the confine- 
ment phase is reflected in the excitation spectrum. In 
order to explore this, we propose two observable proper- 
ties, one optical and another transport. 

Model -- We study nanotubes within the effective- 
mass formalism. There are infinitely many bands cor- 
responding to different modes along the circumference of 
the tube, which we label with n = 0,±1, . . .. The way 
in which the tube is wound is characterized by an index 
v = 0, ±1, which in turn specifies whether a discrete set of 
momenta along the tube circumference intersect the two 
Dirac cones at a =K, K' points in the graphene Brillouin 
zone (y = 0; the "(semi)metallic" case), or not (y = ±1; 
the "semiconducting"), for the half- filled band. In addi- 
tion there is the spin degeneracy for a =f, 4- with the Zee- 
man effect neglected here. Each band is characterized by 
amass hvF^ a {n), where k K ' K (n) = j L (n±<p— v/3) and 
vf the Fermi velocity in graphene fl4j|! We take Hvf = 1 
to be the unit of energy. Here we have introduced a 
magnetic field whose flux passing through the tube is ip 
in units of the flux quantum ipo = ch/eo, which acts to 
shift the discrete set of momenta. Assuming ID electro- 
magnetic fields, the system is described by the extended 
massive Schwinger model with a Lagrangian 
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where F^ = d^Ay — d^A^ is the electromagnetic field 
tensor, tp = (iPr,iPl) T the fermion field with iJ)r,^>l 
the left and right moving components, e = e Q /y/e^ the 
screened charge with a dielectric constant e r , and we 



use the relativistic convention g^ v = diag(l,— 1) with 

7° = (So). 7 1 = (?V) and 4 = A m7 ". The model 
has a spin SU(2) symmetry, which is enlarged to SU(4) 
if K and K' are degenerate. However, perturbations 
can break this into SU(2)xSU(2). When the electric 
fluxes arc confined in ID, the inter-electron potential is 
V(x) = e 2 \x\/2, which enforces confinement of charges. 
The coupling strength is e 2 J{kvf) = 0.709/e r deter- 
mined by the screening factor e r (see e.g. |15l|)- 

Exciton spectrum — Excitons are of central impor- 
tance in understanding the optical properties of carbon- 
nanotubes [la - l20j . Specifically, Kishida et al. have ob- 
served not only the bright excitons, but also the dark 
excitons [ljj, while Wang et al. observed excitons in 
metallic nanotubcs|20|. In the confinement phase, the 
optical spectrum is dominated by excitons, and no two- 
particle continuum exists. Since there is a possibility 
that the lightest fermion mass vanishes, the conventional 
weak-coupling approaches cannot be used, and strong- 
coupling methods are required. Two powerful meth- 
ods are applicable; One is the light-front quantization 
scheme |9l-Hlj. and the other is bosonization [6 8]. For 
studying the crossover from metallic to semiconducting 
tubes, we employ the former, since bosonization for a 
massive model is still an open problemQ . With the light- 
cone coordinates x* 1 = (x + ,x _ ) = (x° + x 1 , x° — x 1 )/\/2 
we can eliminate the dynamical variables except for 
the right-fermionic components ipm (with a shorthand 
i = (n, er, a)) by means of the equation of motion. The 
Lagrangian reads 

L = / dx~C = iy/2 / dx~ ^ : 4>\ R d+ip iR : 



+_ J dx-dy~j+(x-)\x- -y \j+(y ) (2) 



with the U(l) current j M = J2i '■ ipil^ipi '■, 1 + — 
(°fi ), j~ = ( ^ ). The index i runs over infinite num- 
ber of modes with the mass term depending on n, which 
contrasts with the standard SU(iV) massive Schwinger 
model where the mass is common to all i's. The free-field 
expansion is ip iR (x-) = ^i J °° ^^ [b t {k+)e- lk+x ~ + 

d\{k + )e lk x ], where 6^ creating electrons and <& holes 
satisfy a canonical commutation, {bi(k + ),b-(l + )} = 

{di(k + ),d]{l+)} = 2Ttk+8 l3 8{k + -l+). The virtue of using 
the light-front formalism is that the groundstatc, which 
is a confinement phase (CP), is described by the Fock 
vacuum |0) C p = |0) Foc k (with 6 i |0) Fo ck = ^|0) F ock = 0). 
Now let us look at the two-particle excitation, i.e., an 
exciton with a wave function 
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where the integral is restricted to k\ + kz = P. From the 
Lorentz invariance, the exciton wave function satisfies 
the Einstein- Schrodinger equation, 2P~P + \ip) — M 2 \ip), 
where P~ is the light-cone Hamiltonian, P + the momen- 
tum operator for the center of mass momentum P with 
P + \ip) = P\ip), and M the excitation energy ("mass") 
of the boundstate. With the light-cone operators, the 
Einstcin-Schrodingcr equation for the wave function is 
given explicitly as 
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with re-scaled momenta k = k/P, k! = k'/P [21(. This 
is an extension of the 't Hooft-Bcrgknoff cauation|9l[22l|. 
here possessing infinite number of modes labeled by i. 
The last term is the anomaly term, which physically cor- 
responds to a virtual process of exciton pair-annihilated 
into a photon and then regenerated as an exciton. The 
process is intimately related to the photon-exciton cou- 
pling, and for bright excitons this term is nonzero. By 
contrast, the term disappears for dark excitons with 
Ylii^i = 0, an d eqn.flU) reduces to the 't Hooft equation 
for planar QCD[22|. For the optical activity the bright 
excitons must satisfy the condition that tp n -fa = ipnia to 
be a SU(2) spin singlet (note d* , creates a hole with 
spin —u), and, when K and K' are degenerate, an addi- 
tional condition VWk = V'no-K' is imposed to make it a 
SU(2) valley singlet. We note that in eqn. (0} the effect 
of vacuum polarization and self-energy corrections due 
to "meso n" p ropagators jT2| are neglected for simplisity. 
We solve [23 the 't Hooft-Bergknoff equation using the 
basis-function method |10.lll|. 

Let us first look at the excitation spectrum (exciton 
energy M) against the magnetic field for a metallic nan- 
otube with v = in Fig. [TJ We immediately notice that 
the system is no longer metallic due to charge confine- 
ment, namely, the spectrum for the bright exciton has a 
gap, and an excitation continuum docs not exist, cither. 
This is in sharp contrast with the case for the conven- 
tional weak-coupling picture with a 1/r potential (inset 
of (a)), where a continuum exists down to zero energy 
at ip = 0. The spectrum has a periodicity with a pe- 
riod tp = 1, which originates from the mass structure of 
the fermion modes. We label each exciton mode with 
(n, I; a), where a =K,K' is the dominant valley character 
near ip = 0, and I — 0, 1, . . . the exciton quantum number 
that labels the bound state in a trapping potential (see 
Fig. [3](c)). Odd-/ states are parity odd and one-photon 
allowed, while even-/ states are only two-photon accessi- 
ble. For v = the (n,l;K) and (—n,l;K') excitons are 
degenerate due to the valley symmetry, so we can omit 
a from the index. 

What can we learn more about the spectrum? Com- 
parison with the meson spectrum in QCD becomes inter- 
esting. In fact, the bright and dark excitons with low- 
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FIG. 1: (Color online) Two-body excitation spectrum against 
magnetic field for bright (a,c) or dark (b,d) excitons for metal- 
lic (a,b; v = 0) or semiconducting (c,d; v — 1) nanotubes with 
L = 10, e 2 /tvvf = 0.2. Solid (dashed) lines correspond to 
states with even (odd) I. Inset in (a) is a schematic spectrum 
for a system with 1/r potential near <p = 0. 
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FIG. 2: (Color online) Two-body excitation spectrum in zero 
magnetic field against circumference L (a,c; with e /ttvf = 
0.2) or against coupling strength (b,d; with L — 10) for metal- 
lic (a,b; v = 0) or semiconducting (c,d; v — 1) nanotubes. 
Solid (dashed) lines correspond to states with even (odd) I. 



est energies (FigdJ have many aspects shared by the r\ 
and 7T (pion) particles, respectively, in QCD[7|. Inter- 
estingly, the r\ and -k boundstates behave differently in 
the strong-coupling limit; As the fermion mass goes to 
zero, the it mass vanishes, while the r\ mass remains fi- 
nite due to a U(l) anomaly. For QCD there is the Gell- 
Mann—Oakes—Renner re/aiion[2J] , which is a relation, 

1/2 

M v oc m ' ark , between the pion and quark masses. When 
translated into the present problem, the relation applies 
to the lightest dark exciton. For example, as shown in 
Fig. Q] (b,d), the energy of the lightest dark exciton, i.e., 
pion in the QCD context, goes to zero at <p = (p = 0, 1 
in metallic nanotubes and at (p — 1/3, 2/3 in semicon- 
ducting ones. However, the exciton energy behaves as 
M ex \ip — ip~\ a with a power a ~ 0.5 and, since the 
fermion mass is proportional to \cp — (p\, we can regard 
this as a manifestation of the G ell-Mann— Oakes— Renner 
relation in nanotubes, which holds even in 1+1D systems 
where chiral symmetry is 'almost' broken [25| . 

The situation is even more interesting for the lightest 
bright excitons (~ ^-particles) , since they remain massive 
even though the fermions are massless. This can be seen 
in the (0,0) and (0,0;^) states in Fig. Q] (a) (c), which 
have nonzero minima at (p = (p. This is due to a U(l) 
anomaly coming from the pair creation-annihilation pro- 
cess, i.e., the last term in eqn. (0|. The physical picture 
is the following. In an exciton, electrons and holes are 



continuously created and annihilated when the fermion 
mass is small, and a cloud of photon is formed. The 
electromagnetic energy of the photon cloud is finite, and 
contributes to the exciton energy. This process is not re- 
stricted to the ID Coulomb potential, and similar effect 
may take place in metallic single-wall nanotubes where 
excitons with finite binding energy were observed [20J. 

In experiments it often happens that the nanotubes 
have random circumference L. In Fig|2ja,c), we plot the 
dependence of the exciton energy on the circumference. 
The mixing between different modes n is small when L 
is small, since the energy difference between modes is 
ex 1/L, and in the limit L — > 0, the system approaches 
to the pure SU(iV) massless Schwinger model [7|. In this 
limit, the lightest bright exciton (?y-particle) mass is given 
by M v = ^/7Ve 2 /7r, which amounts to A/,, = 0.8944 in 
metallic (y = 0, SU(4)) nanotubes and M v = 0.6325 in 
semiconducting (y = 1, SU(2)) ones for e 2 /ttvf = 0.2. 
The mixing of modes lowers the exciton energy as shown 
in Fig|2] (a,c). This is because the mixing reduces the 
anomaly term, where the (0, 0) mode converges to the 
dark mode in the SU(oo) massless Schwinger model. In 
Figj2] (b,d), we plot the dependence of the spectrum on 
the screened interaction parameter e 2 . Starting from 2k 
(weak-coupling limit), the exciton energy increases as the 
interaction becomes stronger, where the increase is larger 



for larger quantum number I. 




0.0 0.2 0.4 0.6 '0.8 L0 0.0 0.2 0.4 0.6 0.8 1.0 

hi/ hi cr hi/ Hi C1 - 

(c) E/E cr =0 E/E cr = 1/2 E/E„=l 




(d) ® 0- 

£=e/2 -ell ell -ell 

FIG. 3: (Color online) J-E characteristics(a) with a logarith- 
mic plot in the inset, and exciton levels(b) for semiconducting 
nanotubes (for £ = 100) in the weak-coupling, nonrelativistic 
limit, (c) Schematic trapping potential V(x) and eigenstates 
for various values of the electric field, (d) Coleman's half 
asymptotic state when the external field is E — E CI . 

Nonlinear transport and half- asymptotic state — Non- 
linear transport in correlated systems is attracting con- 
siderable interests (see e.g. |26[), and in nanotubes it 
serves as a powerful probe in experimentally detecting 
new states of matter [2| ■ Here let us show that the semi- 
conducting nanotubes with ID interaction should exhibit 
a power-law J-E characteristics in external electric hclds, 
as in the 2D organic systems [l3l| . Unlike the 2D case, 
however, the carriers in the field-induced metallic state 
are only "half- free" . We explain how this has to do with 
the half-asymptotic state predicted by ColemanQ- 

Let us look at a semiconducting nanotube with a finite 
fermion mass k in a strong electric field. Since the strong- 
coupling analysis in finite electric fields is still an open 
issue, here we focus on a weak-coupling approach where 
the U(l) anomaly does not contribute with the bright and 
dark excitons becoming nearly degenerate. In the nonrel- 



ativistic approximation, Dirac particles becomes nonrel- 
ativistic fermions with quadratic kinetic terms (p /2M), 
and the exciton binding problem in an electric field E is 
reduced to solving a ID Schrodinger equation (c.f. eqn. 
(4.1) in 0), 

[-dl + V{x)]i^\x) = (Mi - 2k)^ 1) (x), (5) 

c 2 



V(x) 



x e 



-M/C 



cqEx, 
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where Mi is the exciton energy, V(x) the trapping po- 
tential, and we have introduced an exponential damping 
factor with a cutoff £. We plot the J-E characteristics 
in Fig. EJa) for various values of temperature, where the 
nonlinear current is seen to behave like 



J = <T exp[-A(E)/(2k B T)]E, 



(7) 



as in ref.[T3| with A(E) the activation energy. The cur- 
rent follows a power-law until the conductivity reaches a 
peak at a critical field E = E c ? = eo/2e r . The power 
of the conductivity is proportional to T _1 as in ref. [13[ 
but there is also a cutoff £ dependence. At the critical 
field, one side of the trapping potential for a test charge 
becomes flat as depicted in Fig. [3{c). Then a contin- 
uum spectrum emerges, as seen in Fig. (3Jb) where we 
plot the eigenenergies of the Hamiltonian (eqn. (JSJ) ob- 
tained by gluing two Airy functions. There, Coleman's 
half- asymptotic state - a configuration with alternating 
charges but possibly random displacements (Fig.[3Jd)) - 
has the lowest energy, where the external field and the 
force from surrounding charges balance with each other. 
Strictly speaking, carriers with opposite charges must 
switch their position in order for the current to flow, and 
this violates the half-asymptotic state condition, which 
may modify the simple relation eqn. ([7]) for the nonlinear 
current. This is out of the scope of the present work, but 
will merit further studies. 
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